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Abstract 

Chiral symmetry is consistently implemented in the two-nucleon problem 
at low-energy through the general effective chiral lagrangian. The potential 
is obtained up to a certain order in chiral perturbation theory both in mo- 
mentum and coordinate space. Results of a fit to scattering phase shifts and 
bound state data are presented, where satisfactory agreement is found for 
laboratory energies up to about 100 MeV. 



1 Introduction 



The problem of deriving the interaction potential between two nucleons continues 
to be one of the most fundamental problems in nuclear physics. Early field theo- 
retical work in this area [1, 2, 3, 4] encountered many difficulties, mostly due to 
the non-renormalizability of meson theory. This was followed by more phenomeno- 
logical approaches which utilized empirical forms for the medium and short-range 
parts of the interaction potential [5, 6]. During the last two decades a compro- 
mise approach has been developed in which meson exchange potentials provide 
the medium and long-range parts of the nucleon-nucleon (NN) potential while the 
short-range dynamics is treated phenomenologically [7, 8, 9]. Although the latter 
approaches have achieved very impressive empirical descriptions of nucleon-nucleon 
bound state (deuteron) and scattering data the connection between the nucleon- 
nucleon interaction and the fundamental, underlying dynamics of the strong in- 
teraction remains unclear. It is for this reason that the nucleon-nucleon problem 
continues to be of fundamental interest. 

It has been argued [10] that Regge phenomenology can be extended to low- 
energy nucleon-nucleon scattering with Regge poles leading to a one-boson-exchange 
(OBE) potential where i) the contributions of meson trajectories (including scalar 
e's) are dominated by the particles with lowest spin which couple to nucleons 
with a gaussian form factor and ii) gaussian potentials arise from the Pomeron 
and tensor trajectories. Such a potential in a nonrelativistic expansion has been 
constructed by the Nijmegen group [7] and fits data very well. However, Regge 
cuts are simply neglected. The Bonn group [8] made a serious attempt to include 
multi-boson exchange in the framework of old-fashioned perturbation theory. In 
addition to the OBE of known mesons they included the following: 27r and tt p 
exchange with both nucleons and A isobars in intermediate states, "correlated" 
two-pion exchange in the form of a a' scalar meson, ttctobe exchange (with (Tqbe 
an approximation to 27r, a' and tt p exchanges), and ttlo exchange. Agreement with 
data is quite good. 

Nevertheless, the justification for such approaches in terms of quantum chro- 
modynamics (QCD) remains mysterious. In particular, it is not clear how to 
consistently deal with the exchange of mesons which have masses of the order of 
the typical inverse hadronic radius set by the QCD scale A^qcd- This has led a 
number of researchers [II] to attempt derivations of nucleon-nucleon scattering 
from quark models (either constituent or bag) formulated in terms of some effec- 
tive degrees of freedom which carry the same quantum numbers as the current 
quarks and gluons. Although such models are not derived from QCD either, they 
usually have only a few parameters, most of which are fixed by fitting one-nucleon 
properties. Generally these models produce adequate short-range interactions [II], 
but the long range potential continues to be formulated in terms of pion exchange. 

It seems natural, therefore, to start a treatment of the nuclear force problem by 
recognizing the unique role played by the pion. Although we are largely ignorant 
of the non-perturbative dynamics of QCD at low energies, we know there exists 
an approximate chiral symmetry which is broken by the vacuum. This symmetry 
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restricts the form of the allowed interactions of pions among themselves and with 
other particles. Consequences of approximate chiral symmetry are i) the small 
mass of the pion relative to the QCD scale, A^qcd^ and its subsequent long-range 
contribution to the NN potential and ii) theorems relating processes involving 
different numbers of pions which yield some predictive power. The pion is indeed 
the most important character, besides the nucleon, in the nuclear physics drama. 

The distinguished status of the pion in determining the NN interaction has, of 
course, been emphasized before, particularly by the Stony Brook and Paris groups 
[9]. The coordinate space potential developed by the latter contains: i) a long range, 
"theoretical" part constructed through unitarity, analyticity and crossing relations 
from TTTT and ttN phase shifts, which includes one, two (continuum plus /?, e) 
and partially three pion exchanges (in the form of u) and ii) short range, purely 
phenomenological spin and isospin dependent parts. Both groups evolved from 
this model-independent but parameter-crowded approach to the other extreme, 
the two-parameter Skyrme model. Semi-quantitative success resulted, except for 
the lack of a central, intermediate range attraction [f2]. (For a review of further 
developments, see [f3].) 

What is fundamentally new in the present approach [f 4, f 5] is the development 
of the NN potential within the framework of the general effective chiral lagrangian. 
By considering the most general lagrangian which involves the pion and the nu- 
cleon, and transforms under chiral symmetry as the QCD lagrangian, we divide 
the NN problem into two parts. The first task concerns QCD and its reformula- 
tion in terms of the relevant, low-energy degrees of freedom. The resulting theory 
must have the form of the general chiral lagrangian (because the latter contains 
all the interactions with the correct symmetry), where the coupling constants are, 
in principle, known functions of fundamental quantities like A^qcd and the quark 
masses. In other words, the dynamics of QCD is buried in the couplings of the ef- 
fective chiral lagrangian. Since different models based on QCD represent different 
attempts to capture the essence of this underlying dynamics, they will generally 
differ in the strengths of the low-energy parameters. The second part of the prob- 
lem is to relate the parameters of the effective chiral lagrangian to the measured, 
low-energy NN scattering data and deuteron properties. 

Clearly, we do not attempt to "solve" QCD here, but instead concentrate on 
the second task described in the preceeding paragraph. We start with the gen- 
eral chiral lagrangian with undetermined coefficients. Because chiral symmetry is 
manifest (contrary to most meson-exchange models — e.g. [7, 8]), our approach is a 
priori compatible both with QCD and with all known low-energy phenomenology, 
including tttt, ttN , and scattering, meson-exchange currents, etc. Our scheme 
is model independent in the sense that we do not adopt either a massive meson 
exchange picture or a particular quark model. When a systematic analysis based 
on a chiral lagrangian is carried out for such processes as ttN scattering, a number 
of the unknown coefficients in our model can be determined independently of the 
NN data. In the meantime we keep these parameters free in the NN data fitting 
procedure. 

We do have to make one assumption, that of naturalness, which requires that 
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the parameters be consistent with naive dimensional analysis. With this one as- 
sumption a perturbative treatment of the nuclear potential can be developed that 
is lacking in other approaches. Here the perturbative expansion is in powers of 
momentum divided by a typical QCD mass scale. Up to a given order of expan- 
sion the effective chiral lagrangian specifies precisely the terms which appear in the 
NN potential. Of course, there is no guarantee that the resulting potential will be 
sufficient to describe the data. If a good overall description of the data results, it 
means that the perturbation expansion was carried out to the order the precision 
of the data requires. If, on the other hand, an important phenomenological ingre- 
dient {e.g. scalar isoscalar attraction) is missing, then this might indicate that a 
certain operator or diagram is more important than naively expected. This in turn 
would be indicative of some characteristic dynamic mechanism, and we would be 
learning something about QCD. 

We emphasize that our aim here is not to obtain better fits to the nucleon- 
nucleon data than the already excellent fits achieved with meson-exchange poten- 
tials. We do intend for this approach to help establish a bridge between QCD 
and nuclear physics and to provide a sound model of the nucleon-nucleon po- 
tential whose off shell structure is fixed and which may be used for calculating 
other nuclear processes. In short, the general chiral lagrangian is a useful way to 
parametrize both our ignorance of QCD and our knowledge of nuclear physics. 

General ingredients and properties of effective chiral lagrangians for nuclear 
physics applications are discussed in Sec. 2; the effective chiral lagrangian expan- 
sion used here is presented in Sec. 3. The two-nucleon potential is derived to a 
certain order in chiral perturbation theory in momentum space in Sec. 4 and in 
Sec. 5 is transformed into coordinate space, using a momentum space gaussian 
cut-off. The special techniques required to calculate NN scattering and bound 
state properties with the present coordinate space potential are discussed in Sec. 6 
and the results of fitting the nucleon-nucleon scattering and bound state data are 
presented in Sec. 7. Conclusions are given in Sec. 8. Finally, many details are 
deferred to the Appendices. An initial report of these results was presented in 
Refs. [14, 15]. 

2 Power Counting 

In this work the low energy NN potential is expanded in powers of momentum 
divided by a QCD mass scale. Typical three-momenta Q exchanged in nuclei can 
be estimated as the inverse of the rms electromagnetic radius {rlf^y^'^ of a light 
nucleus. For example, for the triton with (j'lhY^'^ — l-75fm we find that Q ~ 
m^, the pion mass. In QCD the coupling becomes strong and is dominated by 
non-perturbative effects below a momentum scale M that is roughly given by a 
typical hadronic mass, ~ I GeV. Whenever we face such a two-scale problem it 
is useful to separate the corresponding physics by considering an effective, low- 
energy theory which involves only the relevant degrees of freedom, all with small 
three-momenta Q. Such theories can be formulated with a lagrangian that is local 
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(in the sense that it involves only operators containing fields at the same spacetime 
point) and shares the symmetries of the underlying theory, in this case QCD. The 
dynamical information for modes with momenta >Af is contained in an infinite 
set of parameters. 

What then are the relevant degrees of freedom in the case of low energy nuclear 
physics? Unlike the situation at high energies where quark and gluon degrees of 
freedom are indirectly manifest in the data (e.^. jets, deep inelastic scattering, 
quarkonium production, etc.), low energy nuclear physics does not reveal this un- 
derlying QCD structure in any obvious way. Therefore the relevant fields for this 
study should represent mesons and baryons. Clearly, the lightest stable particles 
in each sector should be included. The pion tt has a mass that is small compared 
to Af, and its pseudo-Goldstone boson nature makes it a fundamental ingredient. 
The nucleon N has a mass mjv which is not small but because protons and neu- 
trons comprise the principle constituents of nuclei they must be included. (The 
explicit appearance of the nucleon mass mjv in the effective theory requires care 
as has been discussed previously [16, 17]). The effects of higher mass meson and 
baryon states will generally be suppressed by the inverse of the meson masses or 
by the inverse of the mass differences between the baryons and the nucleon. We 
retain only those mass states for which this factor is much larger than ~ 1/Af. In 
the meson sector, this implies that we do not explicitly keep the /?, cu, etc. , whose 
masses are > 5.5m^ which are closer to M than to m^. In the baryon sector we 
retain only the A isobar which has a mass mA ~ mjv + 2m^, but do not include 
the N* with mass mjv* ~ fnj^ ■\- 3.5m^ nor any other higher mass baryon state. 
The contributions of these additional fields could be included in a similar way as 
is done for the A. The other octet pseudo-Goldstone bosons and the hyperons 
are also omitted. For simplicity we consider only SU{2) X SU{2)^ however our 
treatment can be readily extended to SU (3) X SU (3) to encompass hypernuclear 
physics. 

The requirement that the low energy lagrangian incorporates the symmetries 
of QCD restricts the form of possible interactions involving tt , N and A, but we 
are still left with an infinite set of interactions i with coupling constants (/i, which 
differ in the number of derivatives or powers of pion mass Ji, fermion fields /i, etc. 
If we knew how to solve QCD at low energies, we could calculate these coupling 
constants directly. Since there is no a priori reason for the couplings in the effective 
chiral lagrangian to be small, no a priori perturbation expansion for the infinite 
set of interactions can be formulated. 

We can proceed only by making an assumption of naturalness which means 
that when a coupling constant gi of mass dimension —Si is expressed as gi = 
giM~^' , the dimensionless coupling constant gi will be of order unity. Of course 
this might not be true for all the couplings and this will become apparent through 
phenomenological data analysis. If a coupling constant is found to be anomalously 
large or small, it may require special treatment at low energies, but this may also 
indicate a particular dynamical or symmetry effect at the level of QCD. 

We now have a natural expansion parameter ~ the contribution of 
any diagram being characterized by the power p of the soft momentum Q. We 
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organize our perturbation expansion by counting powers of Q in the same way 
that is done to get the superficial degree of divergence of a graph, where special 
care is taken with baryons due to explicit factors which contain their large masses. 
In the present effective theory it is assumed that all three-momenta Q <^ ttin; 
nucleons and A's are therefore nonrelativistic ^. 

The first task is to organize the expansion in such a way as to eliminate time- 
derivatives of the fermions in interaction terms, since they would contribute large 
factors. This has been done by redefining the fermion fields in terms of velocity 
eigenstates [19], but also more simply by directly replacing the time-derivatives of 
fermion fields using the equations of motion for the fermions [16, 17]. In so doing 
we generate interaction terms that have already been accounted for, which simply 
result in a redefinition of existing coefficients. 

The second task is to distinguish between so-called reducible and irreducible 
diagrams. Reducible diagrams are those which can be separated into two parts 
by cutting through an intermediate state which contains only the initial or final 
particles. This type of intermediate state produces infrared divergences in the limit 
when the baryon kinetic energy is ignored; when it is not, a small recoil energy 
denominator results which makes the overall diagram bigger than expected by a 
factor ^ I. The contributions of these reducible diagrams are automatically 
included by solving the Lippmann-Schwinger or Schrodinger (in the nonrelativistic 
limit) equations of motion. 

The simplest way to isolate these two types of diagrams is to work in the frame- 
work of old-fashioned, time-ordered perturbation theory. Irreducible diagrams are 
those that contain only intermediate states with energies that differ from the ini- 
tial energy by an amount 0{Q). For an irreducible diagram with Vi vertices of 
type z, L loops, C separately connected pieces and Ej = 2 A external fermion lines, 
the power of Q can be conveniently written as 



is called the index of vertex i. Any reducible diagram can be constructed from 

^ Since we do not know a prion what the scale M is exactly, it is not clear how relativistic 
corrections (which are suppressed by l/mjv) compare to 1/M corrections. A rough idea of 
their relative importance can be obtained from the following naive dimensional argument. The 
nucleon-nucleon potential in momentum space can be written as V{p,p') = al{p,p') where I{p,p') 
is some dimensionless function of the initial and final cm. momenta p and p', respectively, 
and a ~ 27r^/M^ if we count powers of 2 and it a la [18]. Substituting this in the Lippmann- 
Schwinger equation we obtain an expansion in aQruN ~ Qitin /M^. A shallow bound-state 
indicates that this series barely diverges, so we estimate that ~ QniN. This estimate is 
admittedly crude and it is not crucial for our approach but it suggests that relativistic corrections 
0(^^) are 0(-^). If M is actually larger, it only indicates that relativistic corrections are 
relatively a little larger than assumed here. 



,y = 4 - A + 2L -2C + J2 A, 




where 



A, = d, + ^-2 
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irreducible diagrams by connecting the latter with intermediate states with energies 
that differ from the initial energy by an amount 0{Q^ /niN) or smaller. 

Here we deal with diagrams involving only two external nucleons. Irreducible 
diagrams are then two-nucleon irreducible; any intermediate state contains at least 
one pion or isobar. The two-nucleon potential is defined as the sum of such irre- 
ducible diagrams, their contributions being ordered by Eq. (I). The full NN scat- 
tering amplitude is evaluated by iterating the nuclear potential in the Lippmann- 
Schwinger equation, or equivalently, by solving (numerically) the corresponding 
Schrodinger equation. 



3 Effective Chiral Lagrangian 

In order to construct a perturbative expansion in Q/M^ Eq. (I) requires > 0, 
for in this case there is a lower bound for p corresponding to diagrams with the 
maximum number of separately connected pieces, no loops and all vertices hav- 
ing Ai = 0. Corrections with higher p are obtained by inserting loops and interac- 
tions with Ai > 0, and decreasing the number of connected pieces. We will show 
that chiral symmetry requires 

A, > 0. (3) 

Here, for simplicity, we work with QCD with only two light fiavors u and d 
with masses m„ and m^, but it is straightforward to include the strange quark. 
In the limit of vanishing quark masses there is an SU{2) X SU{2) ~ S'(9(4) sym- 
metry which is spontaneously broken to SU{2) ~ S0{3). As a result, there exist 
Goldstone bosons whose fields live in the three-sphere ~ 5'(9(4)/S'(9(3) , with a 
diameter that turns out to be the pion decay constant ~ I90MeV. EoUowing 
Weinberg [16, 17] we use stereographic coordinates tt; the covariant derivative is 
then 

The baryons considered here provide the I /2 and 3/2 representations of the spin 
and isospin SU{2) groups. A nucleon N (isobar A) is described by a Pauli spinor 
(a 4-component spinor) in both spin and isospin spaces, the respective generators 
being denoted by |ct(|(t("^/^)) and t[t^^^'^'^). There are also, of course, 2x4 transition 

operators and T, satisfying 

StSf = -(2% - le^jko-k) (5) 

TaT;^ = ]:{Sab - i£abctc), (6) 
D 

which allow us to couple N and A in bilinear terms with spin and isospin transfer I, 
respectively. 
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The effective chirai iagrangian is constructed out of tlie fieids Z)^, N and A 
and tlieir covariant derivatives, 

V^D, = d^D, + lE^ X (7) 
V,N = {d, + t-E,)N (8) 
V,A = (a, + i(3/2) . ^^)A, (9) 



wliere 



E, = ^7rx D,. (fO) 



Tliis is done by considering aii possibie isoscaiar terms and imposing tlie discrete 
spacetime symmetries of QCD, parity and time-reversal. 

Tliat is not all though, because the quark masses break S'(9(4) explicitly. The 
symmetry breaking terms can be written as a linear combination of the fourth 
component of a chirai four- vector and the third component of another four- vector, 
with coefficients |(m„ -|- m^) and |(m„ — m^), respectively. We account for this 
explicit symmetry breaking by including in the chirai Iagrangian all the terms 
constructed out of tt, and A that transform under S'(9(4) in the same way. 
Their coefficients will then be proportional to powers of these combinations of 
quark masses. That is the way the pion mass arises, m\ oc (m„ -|- m^), so each 
power of m„ -|- will count as . For simplicity we neglect isospin breaking 
terms proportional to (m„ — m^). When the latter are included along similar lines 
we begin to understand why isospin violating effects are so feeble in most nuclear 
phenomena [20]. Appendix A presents further details regarding the transformation 
properties of the field representation used here. 

By writing operators that are chirai invariant or that break chirai invariance 
proportional to the quark mass term, we immediately see that all interaction terms 
have Ai > 0; operators involving only pions have at least two derivatives or two 
powers of and nucleon bilinears have at least one derivative. Chirai symmetry 
therefore guarantees a natural perturbative low-energy theory. 

The index of interaction A^ provides a useful ordering scheme for the chirai 
Iagrangian. Below we denote by C^'"'\ referred to as the n-th order Iagrangian, 
the collection of terms with indices A^ = n. We explicitly show only those terms 
relevant for our application. Since we evaluate diagrams only up to one-loop, 
interaction operators with more pion fields or isobars than those exhibited below 
do not contribute to this potential, although they are there in general, in many 
cases to assure chirai invariance. Note also that we eliminate some redundant 
terms by integrating by parts, by using the equations of motion (e.g. to eliminate 
nucleon time-derivatives), and by applying Fierz reordering [21]. 

The lowest order Iagrangian is 



+N[tdo - 2D-^F-H • (tt X iV) - mN]N 
-2D-^F;^gAN{t ■ a ■ V7r)A^ 
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--CsNNNN - -CxNaN ■ NaN 
2 2 

+A[ido - 2D-^F;H^^I^^ -{izxiz)- mA]A 

-2D-^F;^hA[NT ■ {S ■ V7r)A + h.c] 

+ ... (11) 

where gA is the axial vector coupling of the nucleon, Ua is the ANtt coupling, Cs 
and Ct are the parameters first introduced by Weinberg [16, 17], and as usual we 
work in units where he = 1. 

In this work we will also employ terms with more derivatives and powers of tt. 
The first-order lagrangian is 



£(1) = -^D-'NN[iV7ry --k'^ 



B2 

B 



TT 

+ ... (12) 

where the iJ^'s are coefficients of order (9(1/Af); in particular, the last interaction 
term proportional to B-^ contributes to a scalar-isoscalar term similar to the a term 
in meson exchange potentials. The second-order lagrangian is 

1 - ^ A - ^ ^ — 

£(2) = iVV^iV- -i[iV(i-a- V7r)V2iV + V^iVrt -a- V7r)iV] 

2mjv iSr 



N — — 

■-^ViV(i • a - Vtt) • ViV 



-C[\{NVNf + (ViViV)^] - C'^^NVN) ■ (ViViV) 

-iC'^lNVN ■ (ViV X gN) + (ViV)iV • (iVa x ViV)] 

-iC'^NNi^JN ■ G X ViV) - zC^(iVaiV) • (ViV x ViV) 
-{C'^SikSji + CgSiiSkj + CgSijSki) X 



[iV(7fca,iViV(7;ajiV + a,iV(7fciVajiV(7;iV] 

-{C[oS,kS,i + C[,6u6k, + C[^6,,8ki)NGkd,N^GiN 



[d,NGkd,N + a,iVafca,iV]iVa,iV 
+ ... (13) 

where the A' and C[ are additional undetermined coefficients of order (9(l/ilf^). 
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Using this expansion for the lagrangian and the rules for diagrams in time- 
ordered perturbation theory, it is straightforward to construct the interaction po- 
tentiaL Because we ehminated time derivatives in all interaction terms but four 
(those that come together with the pion and fermion kinetic terms in C^^\ and 
the Bi term in -C^^^), and because each of these four terms involves at least two pion 
fields, the interaction hamiltonian is just ( — 1) times the interaction lagrangian, up 
to interactions with more pion fields that do not contribute to the order we are 
working. 

4 The Two-Nucleon Potential in Momentum Space 

We are now in position to calculate any process involving soft pions and non- 
relativistic nucleons. Equations (1), (2) and (3) guarantee that the dominant 
contributions to such processes come from tree graphs with the maximum number 
of connected pieces and constructed out of the lagrangian C^^\ When applied to 
processes with at most one nucleon, this is equivalent to that given by current al- 
gebra. For example, the Weinberg [22] pion-pion and Tomozawa- Weinberg [22, 23] 
pion-nucleon s-wave scattering lengths are readily obtained. But in the late 1970s 
Weinberg [24] pointed out that chiral lagrangians, in addition, provide a framework 
for evaluating corrections to the dominant contributions. The systematic treatment 
of chiral perturbation theory in the mesonic sector began with the work of Gasser 
and Leutwyler [25] and has been extensively studied in the case of SU{'i) X S'[/(3), 
up to Z = 1 and = 2, and including electroweak effects (for an introduction, see 
Ref. [26]). A systematic study of the SU (2) X SU (2) chiral lagrangian for processes 
involving one nucleon was started by Gasser, Sainio and Svarc [27] and is continu- 
ing with the work of Bernard, Kaiser and Meifiner, and many others (for a review 
see Ref. [28]). In principle, the coefficients /i^, Bi and A' can be determined 
from analyses of one nucleon processes once all contributions through one loop are 
evaluated. Unfortunately, this has not yet been done. In Sec. 7 we obtain values 
for all the parameters by fitting low energy nucleon-nucleon data. It should be kept 
in mind, however, that the number of parameters in the present potential could 
be reduced when sufficient information from the one-nucleon sector is gathered. 
For the many nucleon system the present theory is consistent with the empirical 
observation that three-(and more-)body forces are smaller than two-body forces. 
Some of the implications of this result are discussed in Ref. [29]. Furthermore, 
meson exchange currents [30], pion scattering [31] and pion photoproduction [32] 
on nuclei have also been studied in the same approach. For the remainder of this 
work we restrict our study to the two-nucleon system. 

For only two nucleons in the initial and final state A = 2 and C = 1; Fq. (1) 
then simplifies to 

// = 2Z + ^KA. . (14) 

i 

As usual we work in the center-of-mass (cm.) system and denote the initial energy 
by 2mjv + initial (final) momentum by p{p') and define q = p — p' and k = 
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+ as the transferred and average momenta, respectively. Subscripts 1 and 2 
on spin and isospin matrices a and t refer to nucleons 1 and 2. 

The leading order potential (with ;/ = 0) is obtained from the graphs in 
Fig.l and interactions given by in Eq. (11). Note that to this order nucle- 
ons are static, so that their energies in intermediate states are simply mjv and 
the A isobar does not contribute. One obtains [16] the well-known static one-pion- 
exchange (OPE) potential supplemented by contact interactions where, 

= -{^-l^)H,.U ^\-Jf'\U Cs + W,-B. . (15) 

The OPE term provides the longest range part of the NN force, and it is well 
established [33] that it accounts for the higher partial waves in nucleon-nucleon 
scattering and the bulk of the properties of the deuteron, such as its quadrupole 
moment. Of course the NN potential has other sizeable components, including 
a spin-orbit force, a strong short-range repulsion and an intermediate range at- 
traction. Clearly, the lowest order result in Eq. (15) does not account for these 
additional components. A test of the present approach is to determine whether 
higher order contributions yield such features. 

Eirst-order corrections in Q/M [v = 1) also come from the graphs of Eig.l, 
but with one vertex from C^^\ However, there are no suitable vertices in Eq. (12) 
for the tree graphs in Eig.l and we conclude that there are no corrections to the 
leading order potential that are smaller by just one power of Q/M^ i.e. 

y(^) = 0. (16) 

This is a direct consequence of parity invariance. Eor the tree graphs, we could only 
add a power of momentum (or subtract one and add an extra power of ml) to 
but this is actually a three-momentum because we eliminated time derivatives. 
This results in an odd number of three-momenta from which parity conserving 
terms cannot be constructed. 

There are, however, many corrections of second-order, where p = 2. This 
includes tree graph contributions from and a number of one-loop diagrams. 

Eirst, we obtain corrections from the tree graphs in Eig.l where one vertex 
comes from the interactions in in Eq. (13) and the nucleons remain static. We 
also obtain tree level corrections where the vertices are from in Eq. (11) but 
where recoil is included in the intermediate state. Order p = 2 tree level correc- 
tions using two factors from C^^^ cannot be formed because, as we noticed above, 
there are no suitable vertices in Eq. (12). The tree graph 0[{Q / M^] correction is 
therefore given by 

2gA_^ ^ (7i • g (72 • g 

E - T^i^P + q^) 
Jq"^ + ml 



^tree 
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• X A;) + Ce^- (Tiq - 



(17) 



where the A^'s and Cj's are combinations (see Appendix B) of the A' 's and C' 's 
of Eq. (13). The exphcit energy dependent term is discussed in Appendix C. 

Second, there are contributions from the one-loop graphs in Fig. 2 with all 
vertex factors coming from C^^\ (Other one loop graphs only contribute to the 
renormalization of parameters in the lagrangian.) Intermediate states include those 
with two nucleons, one nucleon and one isobar, and two isobars. Denoting 



L0± 

A 

straightforward calculation gives 



{q± ly + 4m2 



(18) 
(19) 



V, 



(2) 



loop,noA 



TT 

+ 4 



1 



(27r)-^ c<;_|_c<;_ \loj^—lo 
d^l 1 



3 



LO^ + CO. 



+ ^^^-^ ] ai- {qx /)(T2 • (? X /) 



(20) 



\^c<;_|_ + uj- to 

for the diagrams of Fig.2a,b,c,d that do not include isobars in intermediate states 

d^l 1 



y(2) 



^^^^ t 



(27r)3 {io+ + io_) {io+ + 2A)(cu_ + 2A) 



18 Fi i-V ^ 1 2j 



J3/ 
(27r)3 
1 



[(?' - /^)' + 2ai • ((fx /)a2 • {qx I)] 
(I I \ 



+ 



loj^lo-{loj^ + c<;_) yc<;_|_(c<;_ + 2 A) c<;_(c<;_|_ + 2 A 



+ 



1 1 



1 



+ 



1 



+ 



1 



2 Ac<;_|_c<;_ \loj^lo- c<;_|_(c<;_ + 2 A) c<;_(c<;_|_ + 2 A 

1 \" 



+ 



(cu+ + 2A)(cu_ + 2A) 



+ (3 - 4ii • t. 



dH 



[(^f_l^Y_2a,{qxl)a2-{qxl)] 
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1 



1 



+ 



1 



+ 2A + 2A)(cu_ + 2A) 

1 1 



+ 



+ 



1 



^0;+ + LO- loj^ + LO- + 2A^ 
1 

+ 



\cu_(cu_+2A) cu+(cu_+2A)y 
for the diagrams of Fig.2b,c,d,e with one intermediate isobar, and 

.(2) 



f21^ 



V, 



loop,twoA 



2h\ 
SIR 



(27r)3 
1 



1 



1 



+ (3 + 2ii • t2 



Lo+Lo- {lo+ + 2A) + 2A) 
dH 



1 1 
CU+ + cu_ 2A 



[{f-f'f + ar-{qxl)a,-{qxl)] 



f27rl 



1 



+ 



loj^lo-{loj^ + c<;_ + 4A) 
CU+ + cu_ + 2A / 1 



1 



1 



0;+ + ijj. 



+ 



+ 2A) + 2A) 

1 



+ 2A)2 + 2A 



f22l 



for the diagrams of Fig. 2c,d,e that have two intermediate A's. 

Finally, we consider corrections of order [{Q jMY] where v = ?>. Again, some 
terms could come from the tree graphs of Fig.l with one vertex from C^^\ but the 
same argument used for V^^^ guarantees that 



K3) 

tree 



(23) 



Other third-order corrections would come from the one-loop graphs of Fig. 2 where 
one vertex is from in Fq. (12). Parity invariance requires the contribution 
from Fig. 2a to vanish, as can be confirmed by explicit calculation, and because 
there are no tzNN couplings in C^^\ the diagrams in Fig.2c,d,e also do not con- 
tribute. Fig. 2b gives 



y(3) 

loop,noA 



4 \f3 



dH 



1_ 13(^-2 _ r^^^_B,{^q^ - /^) + imlB,] 
+ 1652<Ti • {q X f)a2 ■ {q X />i • is} (24) 



for no A in the intermediate state, and 



V, 



(3) 



loop^oneA 



l/hA 

9 



dH 



(27r)3 io+io_ {io+ + io_) {io+ + 2A)(cu_ + 2A) 
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+ 452<Ti • [q X f)a2 ■ [q X l)ti ■ 

+ 6BiAio+io_{q^ - 1^)} (25) 

when there is one. 

Further corrections are of higher order [p > 4). They include i) two-loop 
graphs, like the ones in Fig. 3, that are numerous and harder to calculate, and ii) 
tree graphs with a vertex from C^'^\ which would bring many new undetermined 
coefficients. We do not attempt to include them here. 

The momentum space form of the potential, first presented in [14], facilitates 
a discussion of its structure and the comparison with other models. As usual 
the longest range part of the potential is given by one pion exchange [Fq. (15)], 
including the dominant, static OFF potential first obtained by Yukawa [1], plus 
corrections [Fq. (17)]. The Ai and A2 terms in Fq. (17) derive from the leading 
corrections to the ttNN vertex that arise in an expansion of its form factor in 
powers of momenta over the form factor parameter. The q^ dependence is usual (see 
for example Ref. [8] where monopole and dipole forms are used), whereas the 
dependence is not as common, however it too has been recently considered {e.g. 
Williamsburg model [34]). The other correction to the static OFF potential is the 
energy dependent term in Fq. (17), which arises from the recoil of the nucleon 
upon pion emission. 

The intermediate range parts of the potential are due to two pion exchange (TPF) 
and are determined by parameters F^, /i^, — niN^ B2 and B3. The con- 
tributions from box and crossed box diagrams (Fig.2c,d,e) are standard. The one 
in Fq. (20) {g\ term) was first considered by Brueckner and Watson [2], while those 
with A's in Fq. (21) {g\h\ terms) and Fq. (22) {h\ terms) are due to Sugawara 
and von Hippel [4]. As a check, our results also agree with the appropriate limit 
of the expressions listed in Ref. [35]. But we would like to emphasize that there 
also exist TPF contributions from the "pair" diagrams of Fig. 2a, b that are less 
common. Those in Fq. (20) and the i?3-term in Fq. (24) have also been suggested 
before by Sugawara and Okubo [3], but with arbitrary coefficients. Here the terms 
in Fq. (20) are fixed by chiral symmetry in terms of gA and F^y while the B^ term 
comes from the ttN cr-term. To the same order, we also have in Fq. (24) two new 
terms (i?i,i?2)- The corresponding terms with A in Fqs. (21) and (25) are also 
new. It is important to emphasize that these contributions from the non-linear 
coupling of the pion to the nucleon are a consequence of chiral symmetry and that 
they are not usually included ^ in meson exchange potentials {e.g. Refs. [7, 8]). 
On the other hand, these terms are the only form of "correlated" pion exchange 

^More recently, there has been some interest in the constraints of chiral symmetry to the TPE 
NN force, but limited to the diagrams corresponding to Eq. (20). For example, in Ref. [36] the 
scalar- isoscalar component of these diagrams has been studied, although a different definition of 
potential is considered; Ref. [37] discussed the relevance of energy dependence in OPEP to the 
definition of these TPE potentials; and in Ref. [38], Eq. (20) was examined for the unphysical 
case oi gA = 1- 
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in our potential. The more traditional s-wave correlated TPE (Fig. 3a) is higher 
order in the formalism discussed here. 

The loop integrals in Eqs. (20), (21), (22), (24) and (25) diverge. Moreover, 
iteration in the Lippmann-Schwinger equation of (even the lowest order terms 
in) the potential produces further infinities. Regularization is therefore necessary, 
and counterterms are required to absorb the dependence on the regulator. The 
contact terms [the Cj's in Eqs. (15) and (17)] perform exactly this function. Once 
renormalized, they contain the effect of exchange of higher energy modes and are 
not constrained by chiral symmetry; i.e. all combinations of spin operators and 
momenta (up to second power) that satisfy parity and time-reversal are included. 
This results in spin-orbit (C5), spin-spin and tensor (Ct, C3, C4, Ce, C7), and spin 
independent central (C5,Ci,C2) forces. In order to compare with other approaches 
it will be convenient to "undo" our previous Eierz reordering [21] and rewrite the 
coefficients Ci as 

C, = cP + cPt,-t,. (26) 



5 The Two-Nucleon Potential in Coordinate - 
Space 

Nucleon-nucleon scattering calculations, including those presented here, very of- 
ten use a coordinate space representation. In order to transform the momentum 
space potential in Eqs. (15) - (26) into coordinate space we first have to specify 
the regularization procedure. The use of dimensional regularization here poses a 
problem that we have not yet succeeded in solving: how to iterate the potential 
to all orders in arbitrary dimension. Instead we use a momentum space cut-off 
A < Af , as has been done in other potential models, because it is conceptually and 
mathematically simpler. The form of the cut-off function and the value assumed 
for A are somewhat arbitrary and presumably not very important (see results in 
Sec. 7); variations in the cut-off are compensated to some extent by a redefinition 
of the free parameters in the theory. Again for simplicity, we follow the Nijmegen 
group [7] and assume a gaussian cut-off function exp( — /^/A^), which regulates the 
loop integrals in the potential. In order to further regulate the loops arising from 
the iteration of the potential, we also cut-off the transferred momentum q using 
the same cut-off function, exp(— (f ^/A^). 

All integrals over q and / can be reduced to simpler expressions involving one 
dimensional integrals that can easily be evaluated numerically. We use the formulas 
and techniques presented in Refs. [35, 39] — see Appendix D for details. Only the 
final form is presented here. 

The tensor, total spin, and relative orbital angular momentum operators are 
defined, as usual, by 

O12 = 3 — CTi • (72 , 
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S = ^(<?i + ^2) , 

L = -ifxV, (27) 

respectively. In terms ol these operators and the Pauli matrices r in isospin space 
the present potential can be expressed in terms of the following 20 operators: 

QP-1,...,20 _ . . . ^^^^ . ^^^^ ^^^^^ . ^ . 

L ■ Sti ■ T2, L^, L^Ti ■ T2, L^ai ■ ^2, L^ai ■ a2Ti ■ T2, 
{L ■ {L ■ S^T 1 • T2, S12L ■ S, S12L ■ St 1 • T2, Si2L^ , 
5'i2^Vi • r2, S^2{L ■ Sf, S^2{L ■ SfT^ ■ T2 . (28) 

The NN potential in coordinate space is written as 

20 Q Q2 

~^ or ov^ 

where 

^.(-' E) = V:ir; E) + V^{r- E)^ + V^\r; E)^^ (30) 

is an energy dependent radial operator determined by the radial functions V^°(r; £"), 
Vp{r; E) and V^{r] E). These sixty functions (some vanish) are listed in Appendix 
E. Each consists of a sum of terms with coefficients determined by the parameters 
of the chiral lagrangian, and each term involves at most one one-dimensional inte- 
gral of the functions from Appendix D. They are smooth at the origin thanks to 
regularization. The energy dependence in the radial functions of Eq. (30) is linear 
(see Appendix C). 

The first eight operators, 0'^~^''"'^ , are standard and are accompanied in most 
potentials by radial functions with no derivatives. In this model they receive 
contributions from pion exchanges and contact terms. The next six operators, 
Q'P='Si,...,iA ^ complete the set used in the phenomenological Urbana vl4 potential [6], 
where, = = ^ for p = 9, . . . , 14. What is characteristic of the structure of 
our potential is the presence of first and second derivative terms for p = 1, ... 8 and 
the presence of the other six operators (9p=15,.--,20_ these additional terms 

arise from the 0{k^) dependence in the A2, C2, C4, C7, and recoil correction terms. 



6 Solution of the Schrodinger Equation 

Having obtained a coordinate space representation of the potential the next step 
is to solve the Schrodinger equation numerically. The procedure is standard, but 
care must be exercised with respect to the derivative terms. 

As usual, basis functions of definite total isospin /, total orbital angular mo- 
mentum Z, total spin S*, and total angular momentum J (and its third component 
m) were used; the relative cm. NN wave function was decomposed into a partial 
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wave sum of products of radial and spin-angle functions. By projecting onto the 
spin-angle basis a set of radial Schrodinger equations results which can be written 
schematically as 



where 



R = 



(31) 



2/ir- 



-L{L + l) + J2vi'Hon-E 



fir Y 



(32) 



/i is the reduced mass, and () denotes a matrix element between spin-angle basis 
functions. Spin singlet and triplet L = J channels are uncoupled, so for these 
states i? is a single radial function. For the tensor coupled triplet states with 
Z = Jib 1, i? has two components and quantities X^^^ and X^^^ become 2x2 

matrices. 

In order to eliminate first derivative terms we define R = K(f> where the auxil- 
iary function K is chosen such that (f> satisfies an equation with no first derivatives. 
This determines a differential equation for K which depends on X^^^ and X^'^\ 
given by 



. L- J - " ' (33) 

or 2 L J 

where Det(X(^)) ^ and asymptotically K ~ . The boundary condition on K 
for triplet channels is fixed by further requiring that the two components of 
linearly independent as r — )• oo which results in lim, 
the Kronecker delta function. Function K at finite r was obtained by Runge-Kutta 
integration of Fq. (33). 

The resulting differential equation for (f> is of the form 



are 



where Sij is 



d^<f>{r) 



A{r;E)^{r) 



(34) 



where A{r; E) depends on X(°), X^^\ X(2) and K. The wave function cf){r) satisfies 
the usual boundary conditions; i.e. cf) vanishes at r = and for large r, cf){r) 
matches to the asymptotic wave functions appropriate for scattering or bound 
states. The S'-matrix or binding energy is obtained from the latter boundary 
condition. The NN S'-matrix is expressed in terms of the usual phase shifts and 
mixing angles as in Fq. (7) of Ref. [40]. Fq. (34) was solved numerically for several 
positive scattering energies and for negative values of E to determine the deuteron 
binding energy and other properties. The calculated phase shifts and deuteron 
properties depend on the undetermined parameters in the lagrangian. The cut-off 
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parameter A was fixed and the remaining parameters of tlie iagrangian were varied 
untii an optimized fit was obtained to recent NN pliase sliifts [41] (witli errors from 
Ref. [42]) and measured deuteron properties [43]. 

7 Fitting Results for Phase Shifts and Deuteron 
Properties 

Tlie 26 parameters of tlie model (qa, fiA, Fw, ^i, ^2, -Bi, i?2, -B3, Cs^\ Ct\ 
C}°\ . . . 6*7°^ Cs^\ Ct\ Ci'\...C^^\ see Appendix E) were varied in order to 
optimize the fit to the isospin (from np) and isospin 1 (pp) phase shifts of Ref. [41] 
at 10, 25, 50 and 100 MeV laboratory kinetic energy. All partial wave channels 
with total angular momentum J < 2 were included in the fits. In addition the 
/ = 0, ^Si Di tensor coupled bound state (deuteron) binding energy, magnetic 
moment and electric quadrupole moment were also used to constrain the fit. The 
phase shifts for the J > 2 partial waves are dominated by the OPE potential at 
these low energies and were not used in the fitting procedure. The masses for 
the pion, nucleon and isobar used were = 140 MeV, mjv = 939 MeV and 
mA = 1232 MeV, respectively. The principle results of this study were obtained 
assuming the cut-off parameter A to be 3.90 fm~^ (equal to the p mass). Sensitivity 
to the cut-off parameter is discussed later in this section. 

The recent Nijmegen [41] phase shift solution was selected for fitting; errors 
were taken from the 1994 Arndt et al. [42] energy dependent phase shift analysis 
(solutions CIO, C25, C50 and ClOO). The relative weighting of the chi-square 
contributions from the deuteron properties (binding energy, magnetic moment and 
electric quadrupole moment) and the scattering phase shifts was adjusted so as to 
achieve a suitable balance. The model was fitted to the phase shift parameters 
rather than directly to the NN scattering data since our goal here is to demonstrate 
the capabilities of the effective chiral Iagrangian approach rather than to attempt 
to generate a phenomenological description of data which competes with other 
meson exchange models [7, 8, 9]. 

A grid search using parameters /i^, Ai, A2, 5i, ^2, ^3, Cf \ Cf\ Cf\ . . . Cf^ 
and fitting the / = phase shifts and deuteron properties was initially conducted 
followed by a similar grid search for parameters Cs'\ Ct \ C}^\ . . . 6*7^^ for the 
1=1 phase shifts using the previously optimized values of /i^, Ai, A2, i?2, 
and B3. The OPE gA and parameters were held fixed throughout the grid 
searches. A full, 26 parameter grid search was not feasible due to computational 
resource limitations. After locating a minimum in the chi-square space via the 
grid searches, the fits were optimized by simultaneously varying all 26 parameters 
using the downhill simplex method of chi-square minimization [44]. 

The best fits to the Nijmegen phase shifts with A = 3.90 fm~^ are shown for 
/ = and 1 in Eigs. 4 and 5, respectively. Except for a few of the channels at 
100 MeV, the fits (solid lines) are in quantitative agreement with the phase shifts 
(data points) where the errors from [42] are shown if larger than the data symbol. 
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The results are essentially the same as shown previously in Ref. [15] but these 
new fits are in significantly better agreement with the 25 and 50 MeV ^Pi and ei 
Nijmegen phases than was obtained in this earlier analysis of the older SP89 phase 
shift solution [42]. The L = singlet and triplet scattering lengths are predicted by 
our model to be -15.6 and 5.40 fm, respectively, in comparison with the measured 
values of -16.4(1.9) fm [45] and 5.396(11) fm [46]. The optimized values obtained 
here for the 26 parameters are given in Table I. 

The predicted phase shifts and mixing angles from our model (solid curves) for 
energies from 100 - 300 MeV are compared with the Nijmegen phase shift solutions 
(data points) in Figs. 6 and 7. For most of the partial wave parameters, except 
^Pi, ei and £2, the model predictions and phase shift solutions are in qualitative 
agreement. Because of the low momentum nature of the model, as expressed in 
the explicit [Q /M) expansion, no effort was made to fit the phase shifts at energies 
above 100 MeV. 

The deuteron properties for the A = 3.90 fm~^ model fit are given in Table II 
in comparison with the measured values from Ref. [43]. Included are the binding 
energy, magnetic moment, electric quadrupole moment, asymptotic J-state to s- 
state wave function ratio, and J-state probability. The s- and J-state radial wave 
functions are also shown in Fig. 8. The negative portion of the J-state at small 
radii is also seen in the deuteron wave function of the Bonn potential [8], although 
both the radial extent and magnitude are larger here. We do not claim that 
the short range, high momentum components of our potential are realistic; no 
quantitative significance should be attached to this short-range part of the wave 
function. The depletion of the J-state at small radii, however, contributes to the 
low J-state probability of ~3% which we obtain. Both the quadrupole moment 
and the asymptotic d/s ratio are about 10% too small. 

Also given in Table II are corrected values for the deuteron parameters corre- 
sponding to the potential model reported previously [15]. In these earlier calcula- 
tions the deuteron wave function was computed incorrectly, resulting in erroneous 
values for the calculated magnetic moment, quadrupole moment and J-state proba- 
bility ^. For the corrected values the magnetic moment increased slightly by 1.4%, 
the quadrupole moment increased by 10% and is closer to the measured value, 
while the predicted J-state probability decreased from 5% to 3%. The scattering 
phase shifts, mixing angles, deuteron binding energy and asymptotic d/s ratio in 
Ref. [15] are not affected. 

It is interesting to study the sensitivity of the calculated phase shifts and 
deuteron parameters to the terms in the potential which are a direct consequence 
of chiral symmetry, corresponding to the diagrams in Figs. 2a and 2b. These in- 
clude the first two terms in Fq. (20), the first term in Fq. (21) and the potentials 
in Fqs. (24) and (25) which depend on parameters B2 and B3. To study this 
sensitivity, calculations for all partial wave channels were made in which each of 
the above terms in the potential was individually set to zero. The first two terms 

^We thank Prof. K. Holinde for suggesting there should be a mistake in our earlier value for 
the quadrupole moment. 
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in Eq. (20) and the first term in Eq. (21) have minor effects on the scattering 
phase shifts and mixing angies, however the chirai symmetry terms in Eq. (20) 
significantiy affect the deuteron properties. The potentiais in Eqs. (24) and (25), 
with the vaiues for the parameters B2 and B3 given in Table I, contribute sub- 
stantially to the scattering predictions and the deuteron. This applies to Eqs. (24) 
and (25) individually and to the (i?i, B2) terms and B3 "cr-term" individually as 
well. 

The NN potential model presented here is, admittedly, complicated. To assist 
the reader we show in Eig. 9 the radial potentials for the ^ So channel corresponding 
to the A = 3.90 fm~^ cut-off and the parameter values in Table I. The radial 
potentials 1^°, and W'^ are defined by taking spin- angle matrix elements of 
the coordinate space potential in Eq. (29) where 

{V) = W%r-E) + W\r-E)— + W^{r-E)— . (35) 

or or^ 

Eor coupled partial wave channels the W functions become 2x2 matrices. The 
values for the ^ Sq potentials ^ and W"^ at 50 MeV incident laboratory kinetic 
energy (the explicit energy dependence is weak) are shown in Eig. 9 by the dashed, 
dash-dot and dotted curves, respectively. The units for , and W"^ are MeV, 
MeV-fm and MeV-fm^, respectively. We also define an effective, local potential, 
Vef f{r; E)^ according to: 

Air; E) = ^^^^^ + 2/iKjj(r; E) - 2fiE , (36) 

where A{r]E) was defined in Eq. (34). The effective, local potential for this case 
is shown in Eig. 9 by the solid curve. If the first and second derivative terms in 
Eq. (35) were set to zero then Veff{r] E) would be identical to W^{r] E). The small 
difference between the solid and dashed curves in Eig. 9 is due to the derivative 
terms. 

Eits to the phase shifts and deuteron properties were also obtained with cut-off 
parameter values of 2.50 fm~^ and 5.00 fm~^. The results for the phase shifts 
and mixing angles for the A = 2.50, 3.90 and 5.00 fm~^ potentials are shown 
in Eigs. 10 and 11 by the dashed, solid and dotted curves, respectively. Using 
the A = 2.50 fm~^ cut-off the ^Pi phase shift and ei mixing angle were better 
described than with the A = 3.90 fm~^ cut-off, however poorer fits to the ^6*0, 
and ^E2 phase shifts were obtained. Improved descriptions of the ^ Pi and 
^P2 phase shifts were achieved with the A = 5.00 fm~^ cut-off value compared to 
the A = 3.90 fm~^ results, however poorer descriptions of the ei and £2 mixing 
angles and the ^Sq and ^D2 phase shifts resulted. The corresponding deuteron 
parameter values for A = 2.50 fm~^ and 5.00 fm~^ are also given in Table II. 
Overall we find qualitatively similar descriptions of the NN scattering results and 
deuteron properties for a wide range of cut-off parameters from 2.5 to 5.0 fm~^ 
(corresponding to a mass range from 0.5 to 1.0 GeV). 
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8 Summary and Conclusions 



We derived a low energy nucleon-nucleon potential, from an effective chiral la- 
grangian for soft pions and nonrelativistic nucleons using a perturbation expansion 
in powers of (Q/M). We expressed the potential both in momentum space and 
in coordinate space, solved the corresponding Schrodinger equation in coordinate 
space, and fitted scattering phase shifts and deuteron properties by varying the 
undetermined parameters of the lagrangian. 

In spirit, our approach is similar to that of the Paris group [9] where information 
on pion dynamics was used to construct the longer range parts of the potential, 
while more complicated dynamics was buried in unconstrained, short range parts. 
The fundamental difference between the approach of the Paris group and that 
of the present work is our use of effective field theory, rather than dispersion 
relations. Use of an effective chiral lagrangian not only ensures that our results 
are consistent with other aspects of pion phenomenology (chiral lagrangians to the 
order we use generally agree with data at the 20% level), but more importantly, 
explicitly incorporates the symmetries of QCD and provides a natural perturbative 
expansion. In this way we, like the Nijmegen group [7, 10], develop a potential 
within a theoretical framework, but unlike Refs. [7, 10] we carry out a controlled 
expansion. Our use of field theory and old-fashioned perturbation theory, on the 
other hand, causes our potential to be similar to a low-energy version of the Bonn 
potential [8]. 

The potential in momentum space shares several features with these and other 
potentials. The short range parts have all the necessary spin and isospin structure. 
The pion exchange terms result in contributions that have been considered before, 
but also result in several new terms related to chiral symmetry. Energy dependence 
(which has implications for few-body forces [29]) arises naturally. 

The potential was transformed into coordinate space using a gaussian cut-off 
function. The 0{k^) dependence in the momentum space potential leads to first 
and second derivative terms in the coordinate space representation. Elimination of 
first derivative terms in the radial Schrodinger equation through use of an auxiliary 
function permitted standard numerical methods to be employed. 

We obtained reasonable, qualitative fits to the deuteron properties together 
with quantitative fits to most of the scattering phase shifts up to 100 MeV in- 
cident nucleon kinetic energy. This shows that our approach accounts for the 
principle features of the nucleon-nucleon potential and that these features can be 
naturally understood from the symmetries of QCD. However, the present work 
also makes clear that it is not practical for potential models derived from effec- 
tive chiral lagrangians to compete with more phenomenological approaches, with 
respect to obtaining quantitative descriptions of NN data over a wide range of en- 
ergies. Extension of the present model to higher energies and further improvement 
in the description of data could only result by including higher orders in chiral 
perturbation theory. 
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A Appendix 



Pions are (pseudo)Goldstone bosons of the spontaneous breaking 
5* (9(4) — )• S0{3). They are associated with the broken generators of S'(9(4) and 
therefore hve in the sphere S0{4:) / S0{3) ~ S^. If we embed it in the euchdean 
space, S0{4:) transformations can be viewed as rotations of in planes. For 
example, SU{2)v of isospin consists of rotations in planes orthogonal to the fourth 
axis, while axial SU{2)a are rotations through planes that contain the fourth axis. 

The sphere can be parametrized in a variety of ways, for example with four 
cartesian coordinates {(p^(fi4 = cr} subject to the constraint, 

a' + ^' = \fI (37) 

It is more convenient, however, to work with three unconstrained coordinates; 
therefore we use stereographic coordinates where 

Under an SU{2)v transformation with parameter e, the tt coordinates rotate ac- 
cording to 

(57r = e X TT, (39) 
but they transform non-linearly under SU (2) a with parameter i as given by 

= (l - ^) I + ■ ■ (40) 

A covariant derivative [see Eq. (4)] can be constructed, which is an isospin 1 object, 

6D, = sx D,, (41) 

which transforms under axial rotations as if under SU{2)v with a field-dependent 
parameter, 

6D, = (e X ^) X D,. (42) 
Fermions also transform linearly under the unbroken subgroup 

6N = le-tN (43) 
(5A = ze-i(3/2)A. (44) 

In this case too, it is simplest to work with fields that realize the whole group 
non-linearly, i.e. that transform under axial transformations as if under isospin 
with the same field-dependent parameter as in Fq. (42). In this case 

8N = i{i x^)-tN (45) 

6 A = i{i X ^) • i(3/2)A . (46) 
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It can be easily verified that the covariant derivatives of the pion, nucleon and 
isobar (Eqs. (7), (8) and (9), respectively) are indeed covariant; that is, they 
transform under SU{2) X SU{2) in the same way the fields Z)^, N and A do (see 



A consequence of this is that an isoscalar constructed out of Z)^, A and their 
covariant derivatives will automatically be invariant under the whole SU (2) X SU (2) 
On the other hand, objects that transform under the full group as tensors involve 
also the tt field itself. For example, an S'(9(4) vector can be constructed according 
to 



where its fourth component gives rise to the pion mass term in Eq. (11). 



B Appendix 

Here we list the relations between the A^'s, Cj's of Eq. (17) and the A''s, and C''s 
of Eq. (13): 



Ai = -{A[--A',) 

Ui — — 2 2 

Cs = — Cg — 2^C[2 + C[^) 

Ce = —{C''^ + Cg + -C[Q — C[i — C[^) 



The origin of the explicit energy dependence of the present nucleon-nucleon poten- 
tial is discussed here. One of the nice features of the chiral lagrangian approach is 
that it allows systematic inclusion of nucleon recoil corrections, i.e. energy depen- 
dent terms, as exemplified by Eq. (17) of Sec. 4. Here we give a somewhat general, 
though brief, description of how these terms arise. The systematic inclusion of 



Eqs. (41)— (46)). 




(47) 



1 



C Appendix 
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recoil corrections has recently been shown to result in cancellations between re- 
ducible and irreducible graphs in the three-nucleon problem [29, 31]. This justifies, 
within this approach, certain approximations often made in nuclear physics. 
The Lippmann-Schwinger equation for this case is given by 

T/ rpE± 

where it refer to outgoing and incoming wave boundary conditions, 

Vab = (^b^V^a) , (49) 



H = Ho + V, (50) 

and the labels A, B and C denote quantum numbers for the free many-nucleon, 
pion, and isobar states The energy parameters in Eq. (48) are the sum of the 
individual energies of these particles. 

As is well known, Eq. (48) can be iterated to give the so-called "old-fashioned" 
perturbation theory, represented by the expansion 

T^B = Vab+J2^ac-^ ^VcB + J2 ^^^'7^ F~\^^'^7f + ' ■ ■ ' 

C \i^B — i^C) CD i-C/fi — -C/Cj [^B — ^d) 

where the it and the ie are omitted to simplify the notation. Notice that Vab in 
Eq. (51) is not energy dependent. 

Since we are interested in describing the low energy, nucleon-nucleon potential 
we choose the external particles to be only nonrelativistic nucleons. As in Sec. 2, 
it is convenient to introduce the effective potential as the sum of the irreducible 
diagrams of the series in Eq. (51). In the two-nucleon case this means diagrams 
where there is at least one pion or one isobar in the intermediate states (see Eigs. 1 
and 2). The complete set of diagrams can now be obtained by iterating this 
effective potential where the internal lines are two-nucleon lines [A — )• a, nucleons 
only): 

= Vff,afAE) + E Veff,a,{E) ^ , Vff,,fAE) + ■■■ (52) 

Notice that Vej j^apiE) does depend, by definition^ on the energy E [Ep in the en- 
ergy denominators). To make contact with the nonrelativistic Schrodinger equation 
we recall that, for n heavy nucleons. 



Ep-E^ = Ev"'w + Pi -Ev'^w + p' 



1=1 1=1 



E — ' E ' 1" small corrections 

i^i 2?7^JV ~t 2"^jv 

Ep — Ea + small corrections . (53) 
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Up to small corrections, which can be systematically accounted for, the effec- 
tive potential depends on = Pj^/(2mjv), the nonrelativistic kinetic energy. 
Clearly, in the infinite nucleon mass limit (static limit) this dependence vanishes 

and it is only in the O ^-^^ corrections to the lowest order term that they appear 
(Eq. (f7)). 



D Appendix 

In order to obtain a potential in coordinate space we take Fourier transforms with 
a gaussian cut-off function with parameter A (see [35, 39] for details). With 



erfc(a;) = 



IT Jx 



die 



denoting the complementary error function, we define and use the following func- 
tions: 



/o(r 



Si(r,A) 
S2(r,A) 
Oi(r,A) 



I 



I 



A ) 



Sirr 

hir.m^) - G2{X,r) 



Ar niTy 



^ Mr 

ertc — 

V 2 A 



47r 



-/2(r, m^) 



^° (r,m^) - 



47r 



Mr) 



mi 
47rA2 



2(m^r) 



47r 
3m3 



-gl A I 



ArV 

I 



I + m^r + -[m^ry | e 



-erfc(-^ + ^ 
' 2 A 



1 + 



I - m^r + -{m^ry 



--erfc ('^ + ^ 
2 A 



Trr 



2^4 



3m£ 



m34(r, m^) - {ml + X^^e-^ ^^Tir, ^ml + X^) 
^ml^'cir, m^) - ^(m^ + X^fl^e-'^<i>l{r, ^m^ + A^) 

m^(^;;.(r, m^) - {ml + A2)5/2e-i^(^;;.(r, + A2) 
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plus the integrals: 

^2(r) = ^[h{r,m^)-Hr{r)] , 
where / is any function of A and A = mA — mjv- 

E Appendix 

Here we give the explicit forms of the 60 radial potential functions V^'(r), p = 
1, . . . , 20; z = 0, 1, 2, which appear in the coordinate space version of the potential 
in Eqs. (29) and (30). To save space the following combinations of functions and 
derivatives of functions are defined: 



A(/) ^ 7(2/" + V) 

D.ij) ^ f (r-V) 



£i(/) ^ /' + -/ 

r 

£2(/) ^ /'--/ 



S{f,g) = f"g" + ^f'g' 



r 

2 



n^g) ^ -(/v + /y') + 4/v 
^(/,5) ^ 4/v - -(/y + /V) 

Q(/,5) = (-V-2/" + 2m^/-/o)5 , 

where / = /(r) and g = g[r) are any of the functions defined in Appendix D, and 
a prime denotes differentiation with respect to r. 
We then have (where he = 1): 

= -^{-35l4^'°)[5(/2,F2)] + 65l[i?i5(/2,/2) + i?3m^(/^)1 

TT 
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8 



AO) 



t(0)^ 



-Cf)£i(/o) 



,(0,0) 



2-7^(0,0)^^/21 



9 



-3A5('H2,'H2) + 45(/2,'H2) + -4'R,i''°^[>S(G'2,G2)] 



1 

A 



'5(/2+'Hi,/2+'Hi: 



A ^2 -'O 



25l'^^°)[T(G'2,F2)] 



3AD^{n2)-iT{l2,n2) - -ACHi) + 4'R,i''°)[A(G'2)] 



+ ^A(/2+^i: 



81' 



'[A(G2)]} 



.(0)^ 



AD^{n2)-iT{i2,n2 



1 

A 



^/.lA^7^i^'°)[A(G2)] - ^/.^52A7^i^'°) 



[A(G2)]} 



+ 



+ 



y 1 487r ' 



2^. 



1 

32 L 



i?4^'°)[E2] + -— -^^^'^^[02] 



(l-O)r 



4mjv 



128 
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(1) 



+ 4^4 +3^6 +Y^^7 



5i 



1 



A2ml , ■ 



\ 1 



8F^ yniN"'' ' ^T^gA I 



t(0)^ 



3AD2{n2) + 2V{i2,n2) 

8 



A" 



-D2 



-h\n^l''\D2iG2)]} 



1 



3Cf)£2(/o) 



-^glh\ [aD2{H2) + 2P(/2, 7^2) - ^i^2(/2 + ^1)) 
-h\A'n'-l^'\D2{G2)] + ^/.^i?2A7^i^'°)[i)2(G2)]} 
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^ 47r 



"^^^^ + 8 



4mjv 



+ 
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32mjv 



32^ ("^^'^^^ 



^0 
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Table I: Effective chiral lagrangian potential model parameters for A 
based on the fit to the Nijmegen phase shifts [41]. 



3.90 fm- 



gA 






1.33 




Ha 






2.03 






(MeV) 




192 




Ai 


(10-6 MeV" 


-2\ 


-1.38 




A2 


(10-6 MeV" 


-2\ 


2.44 




5i 


(10-2 MeV" 


-1\ 


0.342 




B2 


(10-2 MeV" 


-1\ 


0.854 




B3 


(10-2 MeV" 


-1\ 


1.77 










I = 


I = 1 


Cs 


(10-4 MeV" 


-2\ 


1.12 


0.135 




(10-4 MeV 


-2\ 


-0.266 


-0.689 


c\ 


(10-9 MeV" 


-4\ 


0.661 


0.381 


C2 


(10-9 MeV" 


-4\ 


3.39 


2.97 




(10-9 MeV" 


-4\ 


-0.330 


-0.0295 


C4 


(10-9 MeV" 


-4\ 


-0.144 


0.453 


Cs 


(10-9 MeV" 


-4\ 


2.10 


-0.910 


Ce 


(10-9 MeV" 


-4\ 


0.281 


0.0998 


Cr 


(10-9 MeV" 


-4\ 


0.581 


1.36 
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Table II: Experimental and effective chiral lagrangian model fitted values for the 
deuteron binding energy [BE), magnetic moment (/id), electric quadrupole mo- 
ment (Qe)} asymptotic d/s ratio and J-state probability (Pd)- 



Deuteron Fit to Nijmegen phase shifts [41] SP89 Fits" 

Quantities Fxperiment^ A = 2.50fm-^ A = 3.90fm-^ A = S-OOfm"^ A = 3.90fm-^ 

BE (MeV) 2.224579(9) 2l5 2^24 2l8 2l8 

fid [i^n] 0.857406(1) 0.863 0.863 0.866 0.863 

Qe (fm2) 0.2859(3) 0.246 0.249 0.237 0.253 

1] 0.0271(4) 0.0229 0.0244 0.0230 0.0239 

Pp (%) 238 2^6 2^0 2.89 



"Corrected values given here for fit in Ref. [15] 
'See Ref. [43] 
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Figure Captions 

Figure (1) : Tree graphs contributing to the two-nucleon potential (sohd 
hues are nucleons, dashed hues pions). 

Figure (2) : One loop graphs contributing to the two-nucleon potential 
(double lines represent nucleons or isobars). Only one time 
ordering is shown for each type of graph. In (d) and (e) we 
only consider those orderings that have at least one pion or 
one isobar in intermediate states. 

Figure (3) : Fxamples of two-loop graphs that are not included in 
our potential. 

Figure (4) : Best fit (solid curves) to the 1 = np phase shifts and 
ei mixing angle from Ref. [41] assuming a cut-off 
parameter A = 3.90 fm~^. Frrors in the phase shifts, 
where shown, are from Ref. [42]. 

Figure (5) : Best fit (solid curves) to the 1 = 1 pp phase shifts and 
€2 mixing angle from Ref. [41] assuming a cut-off 
parameter A = 3.90 fm~^. Frrors in the phase shifts, 
where shown, are from Ref. [42]. 

Figure (6) : Predictions (solid curves) using the A = 3.90 fm~^ cut-off 

and the parameters in Table I in comparison with the Nijmegen 
phase shift solution [41] for the 1 = np phase shifts 
and ei mixing angle to 300 MeV. Frrors in the phase 
shifts, where shown, are from Ref. [42]. 

Figure (7) : Predictions (solid curves) using the A = 3.90 fm~^ cut-off 

and the parameters in Table I in comparison with the Nijmegen 
phase shift solution [41] for the 1 = 1 pp phase shifts 
and €2 mixing angle to 300 MeV. Frrors in the phase 
shifts, where shown, are from Ref. [42]. 

Figure (8) : Deuteron s-state (upper curve) and J-state (lower curve) 

radial wave functions from the present NN potential using the 
A = 3.90 fm~^ cut-off and the parameters in Table L 

Figure (9) : Radial potentials for the ^ So partial wave state at 50 MeV 
using the A = 3.90 fm~^ cut-off and the parameters in 
Table L The potentials W°,W\ and K// 
defined in Fqs. (35) and (36) are indicated by the dashed, 
dash-dot, dotted and solid curves, respectively. The dash-dot 
(dotted) curve corresponds to l^^/fm [W^ /fni^). 

Figure (10): Best fits to the 1 = np phase shifts and ei mixing angle 

from Ref. [41] assuming A = 2.50 fm~^ (dashed curves), 3.90 fm' 
(solid curves), and 5.00 fm~^ (dotted curves). The solid 
curves here and in Fig. 4 are identical. 

Figure (11): Best fits to the 1 = 1 pp phase shifts and €2 mixing angle 

from Ref. [41] assuming A = 2.50 fm~^ (dashed curves), 3.90 fm' 
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(solid curves), and 5.00 fm ^ (dotted curves). The solid 
curves here and in Fig. 5 are identical. 
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